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1. Introduction 

The one-dimensional diffusion in an external potential V(x) is described by the Fokker- 
Planck equation [1] 

J^OM) = J^(M) - 2^-[f(x)P(x,t)], (i.i) 
where f(x) is related to the potential V(x) by 

'M " -55 v «- (1 - 2) 

With P(x,t) = e~ k2t (p(x), equation (1.1) reduces to the time-independent form 

' i2 r 0(x) + 2A[ /(a;)0(a;) ] =A; 2 0(;r) . (L3) 



<ix 2 <ix 

In addition to being a fundamental equation for nonequilibrium phenomena, the Fokker- 
Planck equation is of particular importance in its relation to the Schrodinger equation. 
Setting (f>(x) = e- v ^/ 2 ip(x), we can transform (1.3) into a steady-state Schrodinger 
equation 

-^0*0 + Vs(xMx) = EiP(x), (1.4) 

where E = k 2 , and 

V s (x) = f 2 (x) + f'(x). (1.5) 

Historically, methods developed for the Schrodinger equation have been applied to the 
study of the Fokker-Planck equation. But this should be the other way around as well. 
Theoretically, the Fokker-Planck equation is often more convenient to deal with than 
the Schrodinger equation. With the Fokker-Planck equation we can carry out a more 
systematic analysis, and the results obtained for the Fokker-Planck equation can be 
applied to quantum-mechanical problems described by the Schrodinger equation. In 
this paper, we study the high-energy (large- \k\) behavior of the Green function for (1.3) 
or, equivalently, (1.4). Although we shall mainly work with the Fokker-Planck equation, 
the results of this paper are directly applicable to the Schrodinger equation, too. 

One-dimensional quantum scattering has a long history. In particular, high-energy 
asymptotic behavior of the Green function and related functions have been studied over 
the decades by both physicists and mathematicians. This classical topic has recently 
attracted attention in connection with inverse problems and the theory of integrable 
systems, and this area of research remains active even to the present day [2]-[17]. 

An essential matter in the study of an asymptotic expansion is the estimation of 
the remainder term. In conventional methods, which are mostly based on an integral 
equation, it is necessary to impose strong conditions on the potential in order to control 
the remainder term. For example, the integrability of Vg is often required. For non- 
integrable Vs, different specific methods need to be used. There has not yet been a 
formalism in which the asymptotic analysis of the Green function can be carried out 
systematically and in a unified way for a variety of potentials that may not necessarily 



High- energy asymptotic expansion of the Green function 



3 



be integrable and that may not even be finite as \x\ — > oo. In this paper, we take a 
totally new approach to this problem, and present a new method for the high-energy 
analysis of the Green function which is applicable to a larger class of potentials. 

Our method is based on the analysis of reflection coefficients. Reflection coefficients 
are fundamental quantities in scattering theory, and they serve as building blocks 
for constructing the Green function. The Green function can be expressed solely in 
terms of reflection coefficients for semi-infinite intervals [18]. High-energy behavior of 
the reflection coefficients was studied in a previous paper [19], and formulas for their 
asymptotic expansion were derived there. In the present paper, we apply these results 
to the Green function, and derive the expansion in powers of the inverse wave number. 
The coefficients of the expansion are expressed in a simple form in terms of a linear 
operator, and the remainder term is expressed in terms of transmission and reflection 
coefficients for finite intervals. The validity of the asymptotic expansion can be studied 
by using this expression for the remainder term. The short-time expansion of the Green 
function can also be obtained by this method. 

We assume that the (Fokker-Planck) potential V(x) is a real function which is finite 
for finite x, and which either converges to a finite limit or diverges to infinity (+oo or 
— oo) as x — > ±oo. The quantity k is taken to be a complex number with Imk > 0. 

In sections 2 and 3, we review the relevant results of previous papers. We derive the 
asymptotic expansion of the Green function in sections 4-6, and discuss the conditions 
for its validity in sections 7 and 8. In section 9, the short-time expansion is studied. 
In section 10, we explain how our method can be applied to the Schrodinger equation. 
Various examples are given in section 11. 



2. Reflection coefficients and the Green function 



We define the time-dependent Green function G F (x,y;t) as the solution of 
d 



d 2 d 
dt dx 2 dx ^ X 



G F (x,y;t) = S(x-y)S(t) 



(2.1) 



with the condition that G F (x, y; t) — for t < 0. Physically, G F (x, y; t) is the probability 
density of finding the diffusing particle at position x at time t, under the condition that 
it was initially at position y. We also define its Fourier transform 



G F (x,y;u) = / e lM G F (x,y;t)dt 



(2.2) 



This G F (x, y; oj) is the Green function for equation (1.3) with k 2 = iu. In the same way, 
we define the retarded Green function for the Schrodinger equation as the solution of 



i— + — - V s (x) 



G s (x,y;t) = 5(x-y)5(t) 



dt dx 2 

satisfying the condition G s (x, y; t) — for t < 0. Its Fourier transform defined by 



G s (x,y;E)= / e iEt G s (x,y;t)dt 
Jo 



(2.3) 



(2.4) 
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V s (x) + E 



G s (x,y;E) = 6(x-y). (2.5) 



is the Green function for the steady-state Schrodinger equation (1.4), satisfying 

~ cP_ 

dx 2 

It is easy to see that Gf and Gs are related by 

G F (x, y; U ) = -e-^- v ^ 2 G s (x, y; E = iu). (2.6) 
For convenience, we define 

G(x,y;k) = 2ikG s (x,y;E = k 2 ) (2.7) 

as a function of complex k, and deal with this G instead of Gf or Gs- (We shall refer to 
G as the Green function, too.) Since G(x,y; k) = G(y,x; k), without loss of generality 
we assume that x > y. 

To introduce our basic expression for the Green function, let us first define the 
transmission and reflection coefficients for finite intervals. We consider an interval 
(xi,x 2 ), and define 

r v( Xl ) (x < Xl ) 

V(x) = V{x) ( Xl <x< x 2 ) , /» = -^V(x). (2.8) 
[ V(x 2 ) (x 2 < x) 

Namely, V(x) is identical to V(x) within (xi,x 2 ) and constant outside this interval. We 
consider equation (1.3) with f(x) replaced by f(x): 

-^0(x) + 2-jL\f(x)<Kx)] = k 2 4>{x). (2.9) 

(In general, the left-hand side of this equation contains delta functions at x — xi and 
x = x 2 coming from the derivative of /.) Since f(x) = for x < x\ and x > x 2 , equation 
(2.9) has two solutions of the form 

f e^ x ^- v ^ 2 r(x 2 , xv, k)e- ik{x ~ x ^ x < x u 
MX) = { e-(^) + R r (x 2 , x i; ky>- * x > x 2 , (2 ' 10a) 

J e ifc(x-xi) _|_ ^( X2) Xl - k)e~ ik{x ~ Xl) x < x 1 , 

MX) = { e-^- v ^ 2 r(x 2 , Xl] k)e ik ^ x > x 2 . 

(The factor e ± ^ v ^ X2 " > ~ v< - xl ' ) ^ 2 in front of r is necessary since 0i and <f> 2 are solutions of the 
Fokker-Planck equation, not the Schrodinger equation.) The transmission coefficient r, 
the right reflection coefficient R r , and the left reflection coefficient Ri for the interval 
(xi,x 2 ) are defied by (2.10). 

We can let x 1 — > —oo in (2.10a) or x 2 — > +oo in (2.10b) to define the reflection 
coefficients for semi-infinte intervals: 

R r (x 2 , — oo; k) = lim R r (x 2 , x\, k), Ri(oo, x\, k) — lim Ri(x 2 ,Xi,k). (2.11) 
Let us define 

R r (x,-oo;k) Rijo^k) 
br[X > k) = 1 + R r (x, -oo; k) ' bl{X > k) ~ 1 + Moo, x; k) ' ^ 
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and 

S(x,k) = S r (x,k) + S l (x,k). (2.13) 
It was shown in [18] that the Green function can be expressed in terms of this S as 



G(x,y;k) 



exp 



rx 

ik(x — y) — ik \ S(z, k) dz 
Jy 



(2.14) 



J[l-S(x,k)][l-S(y,k)} 
The logarithm of (2.14) reads 

logG(x, y; k) = ik{x - y) - ik f S(z, k)dz- 1 - log[l - S(x, k)} - l - log[l - S(y, k)}. 

(2.15) 

We shall carry out the analysis of the Green function on the basis of this expression. 



3. High-energy expansion formula for the reflection coefficients 

In this section, we review the formulas derived in [19] for the high-energy expansion of 
the reflection coefficients. Here we deal only with R r . (Corresponding formulas for Ri 
can be obtained in the same way.) First, we define the generalized transmission and 
reflection coefficients, with an additional argument £, as 

R r (x,y; k) - £ 



Rr(x,y;£; k) 



1 - £R r (x,y; k)' 



tr 2 {x,y;£;k) 

l - Z,Rr(x,y,t,; k,y 

VI - £ 2 r(x,y; k) 



Ri{x,V\£\k) = Ri(x,y;£;k) + 
T(x,y;£;k) = 



1 - £R r (x,y; k) 

We also define the operator M which acts on functions of x and £ as 

m 9 (x, o = -f(x) [(e - r 1 ) g{x, o + r'gix, o)] + r 1 J* ^ ^ g ( x , o, 

where g(x,t;) is an arbitrary function, and f\x) is the function in (1.1)— (1.3). 
It was shown in [19] that, for an arbitrary nonnegative integer N, 

R r (x, — oo; £; k) 



where 



c n (x,t) = -(l-?)M nr - 1 f(x), 
1 



5 N (x,£; k) 



K n (x,0 



(2ik) 



N 



t (x, z; f ; k)K N (z, R t (x, z; f ; k)) dz, 



i + 4 



>,0- 



(3.1a) 
(3.16) 
(3.1c) 

(3.2) 



(3.3) 

(3.4) 
(3.5) 

(3.6) 
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In (3.5), K N (z, Ri(x, z; £; k)) is the quantity obtained from K N (z, £) by the substitution 
£ -R;(x, z;(,;k). It is convenient to define 

c n (x, f ) = 1 ^ Cn(x, 0, (3-7) 

so that (3.4) and (3.6) read 

c n {x,t) = -M n - 1 f{x), (3.8) 

K n (x,0 = -{l+^Ycn + i(x,0- (3-9) 
It is easy to calculate c n from (3.8) by using definition (3.2). We have 

ci = -/, C2 = -r + g 3 = -r +f+ vn - fe, 

Ca = -/"' + 5/ 2 /' - (2/ 4 - f - 2ff")i - 3f 2 f'e + f 4 e, etc. (3.10) 

Obviously c n is an (n — l)th degree polynomial in £, whose coefficients consist of the 
powers of / and its derivatives. The K n 's are obtained from (3.9) and (3.10) as 

K = f, K 1 = f- 2f% K 2 = f" -f 3 - Aff'i + 3/ 3 e 2 , 

K 3 = f" - 5ff + 2(2/ 4 - f - 2ff")i + 9ff'e - 4f 4 e, etc. (3.11) 

Equations (3.3)-(3.6) hold for any f(x) as long as c n and 5n make sense. However, 
(3.3) is meaningful as a high-energy expression only if 

lim k N 5 N (x,^,k) = 0. (3.12) 

Using (3.5) with (3.9), we can study the conditions for (3.12) to hold. For simplicity, let 
us assume that f(x) and all its derivatives are monotone for sufficiently large \x\. (This 
condition is unnecessarily strong and can be relaxed, but we make this assumption in 
order to simplify the presentation.) Then it can be shownf that: 

(i) When \k\ — > oo with fixed arg/c in the range < arg/c < tt, equation (3.12) holds if 
(/ is (N — 1) times differentiable and) J^ -1 ) is continuous and piecewise smooth§. 

(ii) When \k\ — > oo with fixed Imk > 0, equation (3.12) holds if is continuous 
and piecewise smooth, and if lim 2 ^_ 00 f(z)e cz = for any positive number c. 

(iii) When \k\ — > oo with Irak = 0, equation (3.12) holds if f( N ~^ is continuous and 
piecewise smooth, and if /(— oo) is finite. 

This means that if the conditions in (i), (ii) or (iii) hold, then (3.12) holds in the region 
(i) e < arg k < n — e, (ii) Imk > e, (iii) Irak > 0, respectively, where e is an arbitrary 
positive number. (We shall always let e stand for a positive quantity.) 

| The proof is given in [19] only for £ = 0, but this is sufficient. (See the comment below equation 
(5.16) of [19].) If (3.12) holds for £ = then it holds for any £. This can be easily shown by substituting 
the expansion of R r (equation (3.3) with £ = 0) into the right-hand side of (3.1a). 
§ This last condition is described as "piecewise differentiable" in [19], but it should properly be 
"piecewise continuously differentiable" or "piecewise smooth" to avoid some pathological situations. 
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4. Expansion of S 

We wish to derive the l//c-expansion of logG from (2.15) by using the formulas given 
in the previous section. To do so, we need to study the expansion of S and log(l — S). 
From (2.12) and (3.1a), we find 

1 



S r (x, k) = lhn i - - — [Rr(x, -oo; f ; k) + f] ■ 



Therefore, the expansion of S r is obtained by substituting (3.3) into (4.1) as 

1 

t ? n I t \ A Sn I t. \ A- • • • A 

2ik 



1 1 

S r (x, k) = —M x ) + 77Ui^s 2 (x) + ■■■ + 7^n s "( x ) + k), 



(2iky 

where, using definition (3.7), 

Sn(2') c n (x, 1), 



(4.1 

(4.2 

(4.3 
(4.4 



a r N (x, k) = lm^ - - 5 N (x, £, k). 
Setting ^ = — 1 in (3.10), we have 

si = -f, s 2 = -f-f\ s 3 = -/" -2//', s 4 = -r-2ff"-f + 2f 2 r + f 
s 5 = -/ (4) - 2ff" - 2f'f" + Aff" + 8/(/') 2 + 8/ 3 /', etc. (4.5 
Substituting (3.5) with (3.1) into (4.4) gives the expression for a r N , 



a r N (x, k) 



where we have defined 



(2ik) 



1 r x 
WJ-c 



t(x, z; k) 
1 + R r (x, z; k) 



K N (z,r]i) dz, 



r 2 (x,z;k) 

rji = Ri{x, z; k) - -r. 

1 + R r (x, z; k) 



(4.6 



(4.7 



The expressions for Si corresponding to (4.2), (4.6) and (4.7) can be derived in the 
same way, using the analogues of the formulas of section 3 for Ri. The result is 



Si(x, k) 



{-2ik) 



tttS^x) + 



1 1 

s 2 (x) + ■■■ + 



\-2ik) 



-2ik) 

t(z, x; k) 
1 + Ri(z, x; k) 



N 



s N (x) +a l N (x,k), 



K N (z,r) r ) dz, 



t 2 (z, X] k) 



(4.8 
(4.9 
(4.10 



1 + Ri(z,x; k) 

The coefficients s n in (4.8) are the same ones as in (4.2). The expansion of S(x,k 
(equation (2.13)) is obtained by adding (4.2) and (4.8): 



S(x,k) 



(2ik) 2 



s 2 (x) + 



(2ik) 



-s 4 (x) + 



+ 



(2ik) 



N 



jSn'(x) 



+ a N (x,k), 



(4.11 



a N (x, k) ee a r N (x, k) + a l N (x, k), 



(4.12 
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where 

, (N (N even) 

N'={ y . 4.13) 

(N-l (N odd) 

Note that (Xjv+i = o~jv for any even number N. 
Obviously, a sufficient condition for 

lim k N a N (x,k) = (4.14) 

|fc|->oo 

is that the following two equations hold: 

lim k N a r N (x, k) — 0, lim k N a l N (x, k) — 0. (4.15) 

|fc|— >oo |fe|— >oo 

From (4.4) we can see that the first equation of (4.15) holds if (3.12) holds. (The limit 
£ — > — 1 does not interfere with the limit \k\ — > oo. We can also show this directly by 
substituting the expansion of R r into the first equation of (2.12).) Therefore, the first 
equation of (4.15) holds under the conditions given in section 3. The conditions for the 
second equation are obvious analogues. Hence it is apparent that (4.14) holds in the 
sector e < arg k < ix — e if f( N ~^ is continuous and piecewise continuously differentiable. 
If, in addition, both /(+oo) and /(— oo) are finite, then (4.14) holds for < arg A; < n. 

5. Expansion of log(l — S) 

From (4.11) we can derive the expansion of log(l — S) as 

-\ log [1 - S(x, k)} = J^j^M x ) + (^jl" 4 ^) + • • • + (2ik) N, 0iN '^ + PN(yX ' ^' 

(5.1) 

where the coefficients 0:4, a^, . . . are expressed in terms of {s n } as 

«2i = ~ — S S 2j 1 S2j 2 • • • s 2jm . (5.2) 
z m=l m {h,-j m } 

(The second sum in (5.2) is over j n = 1,2,3,... for each n (1 < n < m) with the 
constraint Y^=\3n — i-) The remainder term p N of (5.1) can be written as 



1*1/ a N ^ m 



2^m \1-A 



1JV 



1 N '/ 2 

logsi .iv) -Y. n ' 21 



2 ov ^ ^i(2<fc) 



2ri 



(5.3) 



where = 2 X)^] 2 s 2«/ {2ik) 2n , so that 5 = Ajy + o"jv (see (4.11)). The quantity in the 
brackets on the right-hand side of (5.3) is 0(l/\k\ N+1 ) as \k\ — y 00. Hence 

Pn = \cj n [1 + o(l)] + 0(l/\k\ N+1 ) (\k\ 00) (5.4) 
provided that o^r vanishes as \k\ —y 00. Note also that 

p N = p N+ i (N even), p N = p N+1 + ( N odd). (5.5) 



a 
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The coefficients a 2 , «4, ... of (5.1), which are given by (5.2), can be expressed in a 
more compact form. From equations (A. 3) of appendix A, we obtain 

~liL iog[i _ k)] = + ik i Sr ( x > k) - Si ( x > (5 - 6) 

Substituting (4.2) and (4.8) into the right-hand side of (5.6) yields 
1 c) 111 

~2m log[1 " s{x > k)] = wt Mx) + Jm Mx) + Jm Mx) + " ' ' (5,7) 

where we have used S\ = —f. Comparing (5.1) with (5.7) we find that (d/dx)a n (x) = 
s n +i(x) for any even number n. Integrating both sides of this equation gives 

n (x) = J s n+1 (z)dz. (5.8) 

The integral on the right-hand side is uniquely determined in such a way that a n includes 
no additional integration constant when expressed in terms of / and its derivatives. For 
example, s 3 = - (f 2 + /')', s 5 = (2/ 4 + Af 2 f - Iff" + J®)', as can be seen from 
(4.5). In this way, s m is a total derivative for odd m > 3. Hence we have 

a 2 = -f 2 -f, a 4 = 2/ 4 + 4/ 2 /'-2//" + /< 3 \ etc. (5.9) 

If /(— oo) = or /(+oo) = 0, we can write the right-hand side of (5.8) as ff^ s n+1 (z) dz 
or — s n+ ±(z) dz. But (5.8) holds even if neither /(— oo) = nor /(+oo) = 0. 
Equations (5.1) and (5.6) also give the relation between the remainder terms 

Q 

q^Pn(x, k) = ik [<7 r N+1 (x, k) - cr l N+1 (x, fc)J . (5.10) 
6. The 1/fc-expansion of log G 

Substituting (4.11), (5.1) and (5.8) into (2.15), we obtain, for any integer N > 0, 
log G(x, y; k) = ik(x - y) + ^ a i( x , v) + ^ik) 2 ^' v) + ^ik) 3 ^' ^ 

+ ■■■+ J^N a x( x i V) + & N {x, y; k), (6.1) 

where 

rx 

a n (x,y) = - s n+1 (z)dz (n odd), (6.2a) 

Jy 

rx ry 

a n (x,y) = / s n+1 (z)dz+ / s n+1 (z) dz (n even), (6.26) 

rx 

A N (x,y,k) — -ik / a N+1 (z,k)dz + p N (x,k) + p N (y,k). (6.3) 
Jy 

We can also write the right-hand side of (6.2b) as a n (x) + a n (y) with a n given by (5.2). 
For (6.1) to be meaningful as a high-energy expansion, Ajy must satisfy 

lim k N A N (x,y;k) = 0. (6.4) 

|fc|->oo 



High- energy asymptotic expansion of the Green function 



10 



Let us study the conditions for (6.4) using expression (6.3). We first remark that (5.4) 
and (5.5) imply 

Pn = l<J N+ i[l + o(l)] + 0{l/\k\ N+l ) {\k\ oo). (6.5) 
Now suppose that lim|fc|_>oo k N+1 a^ + i = 0. Then, 

lim k N+1 r a N+1 (z, k) dz = 0, lim k N p N = 0, (6.6) 

|fc|— tOC Jy >CO 

and so (6.4) holds. (Since \ctn+i(z, k)\ is uniformly bounded in the interval y < z < x, 
the limit \k\ — > oo and the integral can be interchanged. The second equation follows 
from (6.5).) We know that lim\j t \^ > . 00 k N+1 aN+i = is satisfied under the conditions 
stated at the end of section 4, with iV replaced by N + 1. Hence it follows that (6.4) 
holds for e < arg k < tc — e if is continuous and piecewise continuously differentiable. 
It holds for < arg A; < re if, in addition, both /(— oo) and /(+oo) are finite. 



7. Validity of (6.4) for discontinuous / 



(N) 



The conditions for (6.4) mentioned at the end of the last section are sufficient conditions, 
not necessary ones. Here we show that f( N ) need not be continuous for (6.4) to hold. 

Suppose that f(z) is (M — 1) times continuously differentiable, f ( - M ~ 1 \z) is 
continuously differentiable except at z — 0, and that f( M ^ has a finite jump at z = 0: 

f (M) (z)=C9(z) + ---. (7.1) 

Here 9(z) is the Heaviside step function, and C is a constant. Then f( M+1 ~> contains a 
delta function as f( M+1 \z) = C8{z) + • • •. It is easy to show that K n (equation (3.9)) 
has the form K n = f^ + • • •, where the remaining terms on the right-hand side do not 
contain derivatives of order n or higher. (See (3.11).) Therefore, 

K M+1 (z,C) = CS(z) + ---. (7.2) 

Substituting this into (4.6) and (4.9) with N — M + 1, and then into (4.12), we have 



a M+1 (z,k) = 



C 



(2ik) 



M+l 



J — i 



t(z, w) 



l 2 



1 + R r (z, w) 



S(w) dw 



C 



(-2ik) 



M+l 



t(w, z) 
1 + Ri(w,z) 



C 



(2ik) 



M+l 



t(z,0) 



1 + Rr(z,0) 



9(z) + (-l) 



S(w) dw + ■ ■ ■ 
t(0,z) 



M 



1 + Ri(0,z) 



e{-z) + 



(7.3) 

The quantity in the curly brackets in the last line of (7.3) does not vanish as \k\ — > oo 
when k is real. (See (A. 2).) But its integral does vanish in this limit even when Im k — 0. 
Indeed, using (A. 5), (A.l) and (A. 2) of appendix A we find 





t(z,0) 




1 + R r (z,0)_ 



9{z) + (-1) 



M 



9(-z) \ dz 
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(< 



2ikx 



J2ik max( 



(7.4) 
(7.5) 
(7.6) 



(».o)) 6 {x) + (-1) M (e~ 2iky - e - 2ifemin (^ )) 

+ o(l/|fc|) (|fc| ->■ oo, < argA; < vr). 

Assuming that the part represented by "• • •" in (7.3) is o{l/\k\ M+l ), we have 

a M+1 (z, k) = 0(l/\k\ M+1 ), f a M+1 (z, k) dz = o(l/\k\ M+1 ), 

Jy 

and, from (6.5), 

p M (z,k) = 0(l/\k\ M+1 ). 

From (6.3), (7.5) and (7.6) it follows that A M = o(l/|A;| M ), so (6.4) holds for N = M 
in spite of the discontinuity of f( M \ 

The same argument holds when /( M ) has two or more finite jumps. In summary, 
the following conclusion can be drawn: 

(i) If / is (N—l) times continuously differentiable and /f^ 1 ) is piecewise continuously 
differentiable, and if is piecewise continuously differentiable, then (6.4) holds 
in the angular region e < argA; < n — e with arbitrary e > 0. (Here f^ may 
possibly have a finite number of finite jumps.) 

(ii) If the conditions of (i) are satisfied, and if lim z ^ +00 f(z)e~ cz = lim 2 ^„ 00 f(z)e cz = 
for any c > 0, then (6.4) holds in the half plane Im k > e with arbitrary e > 0. 

(iii) If the conditions of (i) are satisfied, and if both /(+oo) and /(— oo) are finite, then 
(6.4) holds in the upper half plane ImA; > including the real axis. 

(As noted below equation (3.12), we are assuming that f(z) and its derivatives are 
monotone for sufficiently large \z\, at both z — > +oo and z — > — oo.) 

The conditions of (i) are still not necessary conditions. Actually, (6.4) may hold 
even when / is not (N — 1) times differentiable, as will be shown in the next section. 

8. Effects of the jump of /( M ) at higher orders 

Here we show that (6.4) may hold even for N > M + 1 when /( M ) has a discontinuity 
as in (7.1). To this end, let us study Am+i, Am+2, • • • for the case of (7.1). Since 
K M+2 = /( M + 2 ) + • • we have 

K M+2 (z,C) = CS'(z) + ---, (8.1) 

where 8' is the derivative of the delta function. As before, we substitute (8.1) into (4.6) 
and (4.9), and then into (4.12). After integrating by parts, we obtain 

n ! & r w-v n,,\ 1 2 

(?M+2(Z, k) = 




0(-z) + 



(8.2) 



w=0 
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Using (A. 4) of appendix A, we can write (8.2) as 



o- M+2 (z, k) = 



where 



C 



(2ik) M+1 
2C 



r(z,0) 



+ 



(2ik) 



M+2 



1 + Rr(z,0) 



6{z) + (-l) 



t(0,z) 



1 + Ri(0,z) 



0{-z) 



(8.3) 



(8.4) 



1 (M even) 
(M odd) 

In (8.3), we have disregarded the terms of order l/k M+2 which have the form 
Q{z, k)/(2ik) M+2 with Q(z, k) such that J* Q(z, k)dz vanishes as \k\ — > oo. Such terms 
are irrelevant to the following discussion. In addition to C8'(z), the right-hand side of 
(8.1) contains another singular term proportional to f(z)S(z), but the contribution from 
this term can be disregarded for the same reason. We assume that the contribution to 
(8.3) from non-singular terms of Km+2 is o{l/\k\ M+2 ) as \k\ — > oo. 
Now we consider the two cases, y < < x and < y < x. 



\M+2 



e 2ikx _|_ (-^\M e -2iky 



+ o(l/|A;| M+2 ). (8.5) 



(a) y < < x 

The integral of (8.3) can be calculated by using (7.4). For the case y < < x we have 

J y X -M + 2(z,k)dz = J ^ 1 

(Note that the contribution from the delta function in (8.3) is canceled.) From (5.4), 
(7.3) and (A.2) (or from (6.5), (8.3) and (A.2)) we obtain 

C 



p M +i(x, k) 



2(2ik) 



M+l 



& 2ikx + / x 



I A: 



M+l 



), 



p M +i(y,k) = (-1) 



M 



C 



-2iky 



2(2ik) 



M+l 



+ o(l/\k 



M+l\ 



(8.6a) 



(8.66) 



On substituting (8.5), (8.6a) and (8.6b) into (6.3) with N — M + 1, the terms of order 
l/k M+1 cancel out, and we have 

o(l/\k\ M+1 ) (|*|->oo). 



[x,y;k) 



•7) 



Namely, (6.4) holds for N = M + 1. 

The effect of the discontinuity of /( M ) on Am +2 , A« +3 etc can be studied in 
the same way. Singular terms involving derivatives of the delta function such as 
CS",CS"', . . . appear in K M+3 , K M+i , . . ., respectively, but they cause no problems. 
Singular contributions coming from these terms cancel out in the expressions for 
Am+2, Am+3 etc, as we have seen above for Am+i- (The explanation is omitted here, 
but the cancellation between ik f o- N+ i dz and p N is guaranteed by relation (5.10).) As 
a result, the derivatives of the delta function do not produce a term of order l/k N (or 
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lower) in A N . In (6.2) and (6.3), the delta function and its derivatives appear only 
within an integral, and these expressions are well defined|| as long as n, N < 2M + 1. 

However, (6.3) is not well defined for N = 2M + 2 since the square of the delta 
function appears in it. We can easily see that K 2 m+3 contains a term proportional to 
(j(M+i))2 = C 2p In addition, terms proportional to f( M + 1 + i ) fW+i-i) ^ = 1,2,...) 
are contained in K 2M+3 . These terms are ill defined, too, since they produce the square 
of the delta function by integration by parts. To study A 2M+2 , we need to go back one 
step, as we cannot directly use (6.3). By induction, it can be shown that 

M+l 

K 2M+ 2 = f {2M+2) -4£ ]T f(M+i)f(M+l-i) + . . . (g _ 8) 

1=1 

The terms explicitly shown in (8.8) are the ones that give rise to the ill-defined terms in 
K 2 m+3- We substitute (8.8) into the expression for A 2M+1 . Then we can study A 2M+2 
by using this expression together with the relation A 2 m+2 = A 2 m+i — a 2 M+2/(2i/c) 2M+2 . 
After some calculation^ we find that A 2M+2 contains a term of order l/k 2M+2 as 



1 (-1) 



M 



C 2 + ---. (8.9) 



2M+2 (2ik) 2M + 2 ' 2 
So (6.4) does not hold for N = 2M + 2. We can see that (8.9) gives correction to the 
coefficient of order 2M + 2 as 

log G = ik(x - y ) + ±.a 1 + -- + {2ik y M+ 2 [ a *M + 2 + (-l) M C 2 /2] +-... (8.10) 

Thus, the correct coefficient of order 2M + 2 is not o 2 m+2 but a 2 M+2 + (— l) M C 2 /2. 

In summary, when f^ has a finite jump of the form of (7.1), and when the 
jump is located between y and x, equation (6.4) holds for iV < 2M + 1 in the sector 
e < arg/c < 7r — e provided that the part represented by "• • •" in (7.1) is sufficiently 
differentiable. If, in addition, /(— oo) and /(+oo) are both finite, then (6.4) holds for 
< 2M + 1 even for real k. In any case, (6.4) does not hold for N > 2M + 2. Although 
logG can be asymptotically expanded in powers of 1/k even beyond the term of order 
l/k 2M+2 , the coefficients a n given by (6.2) are not correct for n > 2M + 2. In particular, 
«2M+2 is shifted by (— l) M C 2 /2 as in (8.10). (See examples 5-8 of section 11.) 

(b) < y < x 

Similarly, for < y < x we have, by using (7.4), 



/ 

Jy 



|| Needless to say, f* g(z)5 {m \z)dz = (-l) m g^(0) for a < < b. In (6.2) and (6.3), the function 
multiplying is always m times continuously differentiable if n < 2M + 1 (for (6.2)) or TV < 2M 
(for (6.3)). An expression like 9(z)g(z)8(z)dz — 6(0)g(0) may appear in A 2 m+i, where 0(0) = 1/2. 

It turns out that only the first term on the right-hand side of (8.8) is relevant. Anomalous contribution 
comes from where 9(0) appears in formal calculation whereas it should really be 9(0 + e) or 9(0 — e). 
For example, ikf* 9(z)g(z)e lkz dz = —g(0) + g(x)e tkx + o(l) (|fc| -4- oo), where the first term on the 
right-hand side is not —^(0)^(0) = —g(0)/2 but —9(0 + e)g(0). In deriving (8.9), we also use the 
l//c-expansion of R r , Ri and r for finite intervals. (See the comments below (A. 2) in appendix A.) In 
the expansion, we need only to keep track of the terms which are linear in / and its derivatives. 
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PM + i(y, k) = 2(2 ^ M+1 + o(l/\k\ M+1 ), (8.12) 

while the expression for p M+ i(x, k) is the same as (8.6a). Substituting (8.11), (8.6a) 
and (8.12) into (6.3), we obtain 

&M + i(x,y,k) = J^im^ ky + o(l/\k\ M+1 ) (|*| ->oo). (8.13) 

If the limit \k\ — > oo is taken with fixed arg* in the sector e < arg* < 7r — e, 
then Am+i = o(l/\k\ M+1 ) since e 2tky falls off exponentially. The same can be said for 
A M+2 , A M+3 , and so on; the contribution to A N coming from the discontinuity of f( M ^ 
is exponentially small at large \k\ as long as e < arg* < it — e. 

When k is real, however, k M+1 A M+ i does not vanish as \k\ — > oo. So, unlike 
case (a), equation (6.4) does not hold for iV > M + 1 when the limit is taken along the 
real axis. (See example 6 of section 11.) 

The case y < x < can be treated in exactly the same way as (b). It is 
straightforward to extend the arguments of (a) and (b) above to the cases where there 
are two or more such discontinuities. 



9. Short-time expansion of the Green function 



The expansion of G is obtained by exponentiating (6.1) as 

1 



Q — e ik{x-y) 

where 



+ 



(2ik) 



N 



b N + A'j 



N 



K(x,y) = En E a h a h 



m=l m - Y,ji=n 



A' 



N = exp 



N 



E 



i (2ik)' 



+ A 



A? 



N 

E 



(2ifc)' 



(9.1) 



(9.2) 



(9.3) 



o(l/\k\ N ) as 



(Here the prime does not denote a derivative.) It is obvious that A' 
|A;| -)• 00 if A N = 0(1/1*1^). 

The time dependent Green function for the Fokker-Planck equation is obtained 
from G(x,y; k) by the inverse Fourier transformation as 

G F (x, y; t) = - e -F<*)-v<y)]/2 / G{X) y . K ) e -**dw, (9.4) 

Z7T J-00 LlK 

where k is defined by 

ImK>0. (9.5) 



k 2 = itu, 



(When necessary, the integral in (9.4) is to be understood as f^+ ie with positive 
infinitesimal e.) By substituting (9.1) into (9.4) and carrying out the integration term 
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by term, we can derive an expansion of Gp in powers of t. The result is 

(x - yf 



e -[V(x)-V(y)]/2 

G F (x, y\ t) = = exp 

Vint 



where 



At 



(l + g ± t + g 2 t 2 + g 3 t 3 + ■■■ + g N t N + A N ) , 

(9.6) 



9n{ ,V) hx (m-\)\(n-m)\ (x - yY" 1 } 

(See appendix B for the derivation. The expression for Ajv is shown there.) When 
x ^ y, the remainder term A N satisfies 

]im A N (xy;t) ^ Q 
t^o t N 

if (6.4) is satisfied for e < arg/c < n — e. (See appendix B.) For (9.8) to hold, it is not 
necessary that (6.4) hold for ImA; = 0. So this short-time expansion is valid even when 
/(±oo) are infinite. (Here we are assuming that t is real.) 

Note that the right-hand side of (9.7) is not infinite at x = y in spite of the 
appearance of the l/(x — y) 2n ~ m . For x = y we have 

g n (x,x)= hn ^ X \ (9.9) 
yny , j 2 n (2n - 1)!! 

as can be directly calculated. The right-hand side of (9.7) approaches (9.9) as y — > x. 
10. Application to the Schrodinger equation 

From (4.5), we may note that s n for n > 2 can be expressed in terms of the Schrodinger 
potential Vjg (equation (1.5)) and its derivatives 

s 2 = -^s, s 3 = -V£, Si = Vi-V^ s 5 = (2Vg - V£)', etc. (10.1) 

This can be confirmed as follows. From (A. 3c) of appendix A and (2.12), we can show 
that S r satisfies the differential equation 
d 

— S r (x, k) = 2ikS r (x, k)[l - S r (x, k)} + f(x)[l - 2S r (x, k)}. (10.2) 
Substituting (4.2) into (10.2), we obtain 

n 

S 'n = S n+1 — X! S j S n+l-j ~ 2f s n + /<W (10.3) 

Hence it follows that the s„'s satisfy the recursion relation 

n— 1 

s n +i = s' n + '%2s j s n +i- j (n>2), (10.4) 

J'=2 

where we have used si = — /. We can obtain s 3 , s 4 , S5, . . . from this recursion relation, 
starting with s 2 = — Vs- Therefore, any s n (n > 2) can indeed be expressed in terms of 
Vs and its derivatives. Substituting (10.1) into (6.2), we have 

ai = T ^s(^) dz, a 2 = -V s (x) - V s (y), a 3 = - [* \v s 2 (z) - V£(z)} dz, 
J y Jy 

a A = 2V s 2 (x)-V^x) + 2Vi(y)-V^y), etc. (10.5) 
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Substituting (10.5) into (6.1), and returning to definition (2.7), we can write 

1 r x 

log G s (x, y;E) = - log 2i - log k + ik(x - y) + — J V s {z) dz 

~Wf WW + - (^3 f [^M - VffW] dz + • • • . (10.6) 

The recursion relation (10.4) is familiar in soliton theory. The quantities s n obtained 
from it are, apart from the sign, identical to the conserved charge densities for the 
KdV equation [14, 20]. It is known that these quantities also appear in the asymptotic 
expansion for Jost solutions [16, 21]. Our results for the expansion of the Green function 
are valid even when Jost solutions do not exist. 

So far, we have been assuming that the Schrodinger equation was derived from the 
Fokker-Planck equation (1.3). Let us now go in the opposite direction. To transform 
a given Schrodinger equation into a Fokker-Planck equation, it is, in general, necessary 
to shift the energy level. Let ipo(x) be a solution of (1.4) with E = Eq. Then the 
Schrodinger equation (1.4) is equivalent to the Fokker-Planck equation 

"^(*) + 2^[/(aM*)] = ( 10 - 7 ) 

where 

0(x)=^ o (x)^(x), /(x) = £logVo(a:), p=^k 2 -E , k 2 = E. (10.8) 

The relation between Vg and / is now 

V s (x)-E = f(x)+f'(x). (10.9) 

So (1.5) is a special case of (10.9) with Eq = 0. If we want f(x) to be real and finite for 
any finite x, we need to choose E such that ipo(x) > for any finite x. We may take 
E to be the ground state energy if Vs has a bound state. 

Equation (10.6) was derived by using (1.5), which corresponds to E = 0. However, 
since the Schrodinger equation (1.4) is invariant under the replacements Vs — > Vs — Eq 
and k — > p, expansion (10.6) is valid even when Eq ^ 0. We can check this by an 
explicit calculation. Suppose that Eq ^ 0. Then it is obvious that a correct expression 
for logGs is obtained by making the replacements Vs — > Vs — E and k — > p in (10.6) as 

log G s (x,y;E) = -log2i - logp + ip(x - y) + — [ [V s (z) - E ]dz 

lip Jy 

' [Vs(x)+V s (y)-2E } + .... (10.10) 



(2ip) 2 

If we rearrange (10.10) into an expansion in powers of 1/k by substituting 

i i ; -^0 E E El I I E Q 

to = logA:-^- iF + ---, p = Jfc-_-_ + ..., _ = _ + _ + ..., 

and so on, then it becomes (10.6). Thus, (10.6) holds for Eq ^ as well. Hence we 
know that (6.1), too, is valid for Eq ^ 0. (See example 8 of the next section. Note 
that G is not 2ipGs but 2ikGs-) Unlike the coefficients a n , the remainder term Ajy of 
(6.1) is not expressed solely in terms of Vs. But the conditions for the validity of (6.4) 
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remain unchanged when E Q ^ 0. We can understand this by writing (10.10) with the 
remainder term, and then turning it into an expansion in powers of \jk as above. 

Is is not difficult to interpret the results of sections 7 and 8 in the language of the 
Schrodinger equation. We see from the second equation of (10.8) that the conditions 
on / can be interpreted as the conditions on ip . For example, /(oo) takes a nonzero 
finite value when the ground state wave function decays like e~ cx as x — > oo. The 
differentiability conditions on / can be directly related to the differentiability of V s . 
Namely, / is n times differentiate if Vs is (n — 1) times differentiate. 



11. Examples 

Let us consider some simple potentials for which the exact Green function can be 
obtained, and compare the exact expression with our results for the expansion 

, ogG (,, !/ ;« ; ) = «(x- y ) + |L + _|_ + _^_ + .... („.!) 

The coefficients a n are obtained from (6.2), (4.3) and (3.8). Explicit forms of the 
quantities s n (equation (4.3)) are shown in (4.5) for n < 5 . (Alternatively, we may 
use expressions (10.5) in terms of Vs.) We omit the derivation of the exact expressions. 

Example 1. V(z) = 2z, f(z) = -1, V s (z) = 1. 

As the simplest example, let us first consider a linear potential. Substituting / = — 1 
into (4.5), and then into (6.2), we obtain the coefficients Clfi clS 

ai=x — y, a2 = — 2, a 3 = — x + y, = 4, etc. (H-2) 

The exact Green function for this potential is 

k 



G(x, y; k) = -j== exp [W k 2 - l(x - y)\ . (11.3) 

It is easy to see that (11.1) with (11.2) is the correct expansion of the logarithm of (11.3) 
for Im/c > 0. Since the exact logG(/c) — ik[x — y) does not have any singularities in 
\k\ > 1, the right-hand side of (11.1) is a convergent infinite series for \k\ > 1. 
The short-time expansion (9.6) with g n calculated from (9.7) and (9.2) reads 

(x-yf 

The series in parentheses on the right-hand side is Y^=q{— t) n /n\ = e _t , so this expansion 
is convergent for any t. 



G F (x,y;t) = e (x y) exp 

\/Ant 



t 2 t 3 \ 



Example 2. V(z) = z 2 , f(z) = -z, V s (z) = z 2 -l. 

Next, we consider a parabolic potential. This example satisfies conditions (i) and (ii) of 
section 7 for any N but not (iii). From (6.2) and (4.5), we obtain 

ai = — (x-y), a 2 = 2 - (x 2 + y 2 ), 

^^_x*-jf + 2(x 3 -y 3 ) +x _ y ^ a4 = 2 ( x * + y *)-4(x 2 + y 2 ), etc. (11.5) 
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MlogG-ifcQc-y)] Ke[\ogG-ik(x-y)] 




Re k Re k 



Figure 1. The imaginary and real parts of log G(x, y; k) — ik(x — y) for the potential 
V(z) — z 2 (example 2), with x — 1 and y = 0, (a) plotted as functions of |fc|, with 
arg/c fixed at 7r/4; (b) plotted as functions of Refc, with ImA: fixed at 0.75. Solid lines: 
the exact values. Dashed lines: the expansion (11.1) to order l/k N (N = 2 and 8). 



The exact Green function for this potential can be expressed as 
G{x, y; k) = _ , ^YT771 ZrsM x > k )^-(v, fc )> 



where 

^±(x,k) =e~ x2/2 



2 r(-f)r(|-f 



r(-¥)^(-¥ ) |;^)T2,r(|-f)F(|-f 



3. 2 

' 2 > 



;ii.6) 



;n.7) 



Here F is the gamma function, and F is the confluent hypergeometric function defined by 

F(a, r, z) = En=o a XX%X+n-i) ^ zn ■ We can check that ( 1L1 ) with ( 1L5 ) is the correct 
asymptotic expansion of the logarithm of (11.6) when arg/c is fixed in < argfc < 7r 
(figure 1(a)). Now (11.1) is divergent as an infinite series. This asymptotic expansion 
is also valid when Imk > is kept fixed as \k\ — > oo (figure 1(b)). When Imfc = 0, 
however, (11.1) is not valid since log G(k) oscillates and does not vanish as \k\ — > oo. 
The short-time expansion (9.6) now reads 



—x 2 + y 2 



(x 



yf 



(l + g x t + g 2 t 2 + g 3 t 3 



G F (x,y;t) = - 7 =exp 

\/47rt 



2 At 

where the g n 's are obtained from (9.7), (9.2) and (11.5) as 

-1/2 2\ ' ° On 1 

gi = i - o(z + xy + y ), fi-2 



:n.8) 



6 



(x 2 + xy + y 2 ) 



It 



t 2 i , 2\2 

[x +xy + y), 
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93 = ~\ + ^xy + ^{x 2 + xy + y 2 ) + -^{x 2 + xy + y 2 ) 2 - ^{x 2 + xy + y 2 f. (11.9) 
The exact time- dependent Green function has the well-known form 

<?f(*,V;*) = ( "I 1 , , . V^exp f_ [^-^xp(-2t)]n 

V ' y ' ; V^l 1 - ex P(- 4 ^)]/ V l-exp(-4t) / V ; 

The expansion of (11.10) indeed has the form of (11.8) with (11.9). This Gp(t) has 
singularities in the complex t plane where exp(— At) = 1. The nearest singularities to 
the origin are at t — ±ni/2. So, the infinite series (11.8) is convergent for t < ir/2. This 
is a typical case where the short-time expansion is convergent for small t although the 
high-energy expansion, from which it was derived, is divergent. 

Example 3. V{z) = e z , f(z) = -e 2 /2, V s {z) = (e 2z /A) - (e 2 /2). 

This exponential potential satisfies conditions (i) of section 7 for any N but not (ii) or 
(iii). For this potential we have 

ai = ~\ if - ey ) + l (f X - **) . °2 = \ (e x + e") - \ (e 2x + e 2 ^) , 



1 , , ... 3 



a 3 = -- { e* - e «0 + | (e 2 * - e 2 *) + 1 (e 3 * - e 3 *) - 1 (e 4 * - e 4 *) , etc. (11.11) 
The exact Green function has the form 

Z7T r 

G(x,y;k) = -- — - x+ (x,k) + e" k x-(x,k) X -(v,*), (H-12) 

2cos(«7r/ej L J 

where x± are defined in terms of the Bessel functions as 

ie x l 2 r n 1 

x± (x,k) = — =- J v± {-ie x /2)+iJ_ v J-ie x /2) , u± = -±ih. (11.13) 



2 



Equation (11.1) with (11.11) gives the correct asymptotic expansion when arg/c is fixed 
in < arg/c < n (figure 2). However, this asymptotic expansion is not correct when 
\k\ — > oo with fixed Irak, irrespective of whether Irak > or Imk = 0. The short-time 
expansion (9.6) with g n calculated from (11.11) is shown in figure 3. 

Example 4. V(z) = 2 log cosh z, f(z) = — tanh^, V s (z) = 1 — 2sech 2 ^. 

This V(x) tends to +oo linearly as x — > ±oo. Both /(+oo) and /(— oo) are finite, and 

conditions (iii) of section 7 are satisfied for any N. In this case, we have 

a\ = x — y — 2(tanhx — tanhy), a 2 = (sechrr) 2 + (sechy) 2 — (tanhrr) 2 — (tanhy) 2 , 

4 8 
a 3 = —x + y + -(tanhrc — tanhy) + -[(sechx) 2 tanhx — (sechy) 2 tanh?/], etc. (11.14) 
o o 

The exact G for this potential is 

C( h\ k r(q)r(|-a)r(/3)r(|-/3)cosa7rsin/37r rj + (x)ri_(y) 

[X,yj ' 27rv / F^T V(a - £)r(/3 - a) sin[(/3 - a)ir] cosh x cosh y ' 1 ' } 
where 

i]±(x, k) = F (a, (3, |; — sinh 2 x^j 

r(i + «)r(i-/3) /ii, 9 \ , 
t2 r(a)r(i -/3) F ( a + ^ + ^' i ; ~ sinh x ) ' ( 1L16 ) 




Figure 2. The imaginary and real parts of log G(x, y; k) — ik(x — y) for the potential 
V(z) = e z (example 3), (a) plotted as functions of |fc|, with arg/c fixed at 7r/4; x = 0.8, 
y = 0; (b) plotted as functions of x, with k = 2.5 exp(i7r/4), y — 0. Solid lines: the 
exact values. Dashed lines: the expansion (11.1) to order l/k N (N = 2 and 8). 




(a) (b) 



Figure 3. The series 1 + X^n=i 9 n ^ n f° r * ne potential V(z) = e z (example 3), plotted 
as a function of t with various N. (a) x — 0.7, y — 0; (b) x — 3, y = 0. 



with a = |(— 1 — i\/k 2 — 1) and /3 = |(— 1 + i\fk 2 — 1). Here F is the hypergeometric 
function defined by F(a, /3,j; z) = r( ^ ( p } (/3) E^° =0 r( "r" 7 +i / ) + " ) Since /( ±0 °) are 

finite, expansion (11.1) with (11.14) is valid even when Irak = (figure 4). As in 
example 1, the infinite series (11.1) is convergent for \k\ > 1. 




Figure 4. The imaginary and real parts of log G(x, y; k) — ik(x — y) for the potential 
V(z) = 2 log cosh z (example 4), (a) plotted as functions of real k, with x = 0.5, y = 0; 
(b) plotted as functions of x, with k — 1.2, y — 0. Solid lines: the exact values. Dashed 
lines: the expansion (11.1) to order l/k N . (Since k is real, N = 2, 4, 8 are the same as 
TV = 1, 3, 7, respectively, for the imaginary part.) 



Example 5. V{z) = 2\z\, f(z) = 1 - 26 (z), V s (z) = 1 - 2S(z). 

This is an example where f(z) has a jump at z — 0. This belongs to the case of (7.1) 

with M = and C = —2. For y < < x, the exact Green function is 



i + VA; 2 - 1 
G(x, /e) = exp 



i\/fc 2 — i(x — y) 



11.17) 



x — u 

y 3 



+ 0(l/|fc| 5 ). (11.18) 



and hence we have, as \k\ — > oo (0 < argA; < 7r), 

log G(x, y; k) = ik(x - y) + —(a: - y - 2) - 7—^ 

2ik (2ik) 6 

On the other hand, a\ and 02 are obtained from (10.5) as 

a x = f X [1 - 2<J(z)] dz = x-y-2, a 2 = -1 + 2<5(x) - 1 + 2<J(y) = -2. (11.19) 

Comparing (11.19) with (11.18), we can see that a x is the correct coefficient of the 
expansion but a 2 is not. As shown in section 8, the coefficient a 2 needs to be corrected 
by C 2 /2. Since a 2 + C 2 /2 = 0, we can see that (8.10) indeed agrees with (11.18). 



Example 6. 

e z 



V(z) 



l + z 




Va(z) 



(e 22 /4) - (e 2 /2) (z < 0) 
1/4 (z > 0) 
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Figure 5. The graphs of \k N A^(x, y; k)\ for examples 6 and 7, as functions of real k. 
(N = 3 and 4.) (a) Example 6; x = 0.5, y = -0.5. (b) Example 7; x = 0.5, y = -0.5. 



This example belongs to the case of (7.1) with M — 1 and C = —1/2. Now f(z) is 
continuous and f'(z) has a jump at z = 0. For y < < x, the exact G has the form 

r( w -2y/2ike iKx X -(y,k) 

lX ' V ' ] ~ (K - «/_) J„_(-i/2) + (ilf + J-„.(-i/2) 1 UJ 



with K = ^Jk 2 — (1/4), where x~ and i/± are defied by (11.13) . As shown in figure 5(a), 
we can check that (6.4) holds for N < 3 but not for N > A. It can be seen that the 
limit of |fc 4 A 4 | as \k\ -> oo is 2 _5 C 2 ~ 0.0078 as predicted by (8.9). The coefficients of 
the expansion for y < < x are obtained from (6.2) as 

a x = i (-e 2y + 46^ + 2a; - 3) , a 2 = -~ (e 2 ^ - 2e^ + l) , 

a 3 = — (3e 4j/ - 16e 3y - 72e 2 ^ + 96e^ - 12x - ll) , 
3 192 V / ' 

a 4 = - ( e 4y - 4e % - 4e 2?/ + Ae v + l) , etc. (11.21) 
8 ^ ' 

Since both /(— oo) and /(+oo) are finite, expansion (11.1) with (11.21) is correct to 

order 1/k 3 for ImA; > 0. The correct expansion to order 1/fc 4 is obtained by adding 

-C 2 /2 = -1/8 to a 4 as in (8.10). (See figure 6(a).) 

For < y < x, the exact Green function is 

G(x, y,k) = ^ [e-^y + ^-e lK ^ e^, (11.22) 

where 

A± ^{ 1±V ~k) J "-H/ 2 ) " { lTl Jc) J -^H/2)- (H-23) 
The a n obtained from (6.2) for < y < x are 

<H = -(x-y), a 2 = --, a 3 = -— (s - y), etc. (11.24) 
4 2 lb 

We can show that the quantity A^/A + is 0(l/|fc| 2 ) as \k\ — > oo (0 < argfc < it). 

Therefore, from (11.22) we can see that (11.1) with (11.24) is correct only up to order 
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Im(logG)-^-y) Re(logG) 




Figure 6. The imaginary and real parts of log G(x, y; k) — ik(x — y) for the potential 
of example 6, plotted as functions of real k; (a) x = 0.5, y = — 1; (b) x = 2, y = 1. 
In (a), the dashed curve labeled " N — 4" shows the expansion to order 1/k 4 with the 
correction term — C 2 /2 added to as in equation (8.10). In (b), the envelope of the 
oscillation of Re(logG) falls off like 1/k 2 . 



1/k when Imk = (figure 6(b)). For e < axgk < n — e, this asymptotic expansion is 
correct to any order since e tKy in (11.22) vanishes faster than any power of 1/k. 



Example 7. V(z) = (sgnz) 



,1/2 



, f(z) = -ui + \z\y 



-1/2 



" 1 + i(sgnz)(l + b| 



-3/2 



V B (z) = &(l + \z\ 

This is another case where f'{z) has a finite jump at z — 0. The exact Green function 
for y < < x can be expressed in terms of confluent hypergeometric functions as 

2ik(+(x,k)(-(y,k) 



G(x,y; k) 



CUo,A0C-(o,A0-Cf(o,A0C.(o,*) 



;n.25) 



where 

(±(x,k) = e ik{ - l±x) {q 1/2 Y{q) 



F (q, |; -2ik(l ± x)) T (q + 1, §; -2ifc(l ± 



± r (g + |) F (g + i, i; -2iAs(l ± x)) =F (g + |, |; -2ifc(l ± a;] 



;n.26) 
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with q = i/(32k), and (±{x, k) = {d/dx)(±(x, k). For y < < x, the coefficients of the 
expansion calculated by our method are 



at 



a 2 



i [log(i + x) + log(l -y)]-- A [(l + x)- 1 ' 2 

lb 4 L 



1 

16 



(l + x)- l + (l-y) 



-i 



(1 + x) 



-3/2 



y)- 1/2 
-y)- 3/2 



etc. 



(11.27) 



This, too, is a case of (7.1) with M = 1, and so (6.4) holds for N < 3 but not for 
iV > 4. Expansion (11.1) is now correct to order 1/k 3 for ImA; > (figure 7). We can 
see from figure 5(b) that |/c 4 A4| indeed approaches the predicted value 2~ 5 C 2 ~ 0.0020 
as \k\ oo. (In this case, C = 1/4.) 

Example 8. Vs(z) = \z\. 

Here we consider the case where Vs is given, and where E ^ (see section 10). For 
y < < x, the coefficients a n are obtained from (10.5) (or (10.4) and (6.2)) as 



Or, 



a 2 



-x + y, 



-(x 4 + y 4 ) -5(x-y), 



a 3 ■■ 
16 



2-^(x 3 -y 3 ), a A = 2(x 2 + y 2 ), 



«6 = -yi^ 3 -y 3 ) 



+ 10, etc. 



The exact Green function for y < < x is 

ikki(x - k 2 )M(-y - k 2 ) 



G(x,y;k) 



;il.28) 



(11.29) 



Ai(-P)Ai'(-A: 2 ) 

where Ai(z) is the Airy function, and Ai'(z) is its derivative. From the second equation of 
(10.8) we have f(z) = (d/dz) log Ai(|z| —E ), where E is the smallest number satisfying 
Ai'{-E ) = 0. (Numerically, E ~ 1.019.) It is easy to see that /"(0+) - /"(0-) = 2, 
so that this is a case of (7.1) with M = 2, C = 2. Expansion (11.1) with (11.28) is 
correct to order 1/k 5 for Imk > e. (Conditions (ii) of section 7 are satisfied since f(z) 
behaves like |z| 1//2 as \z\ — > oo. This expansion is not valid for Imk = 0.) We can also 
check that the correct coefficient of order 1/k 6 is not a 6 but a§ + C 2 /2 = a$ + 2. 
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12. Summary and remarks 

In this paper, we studied the high-energy asymptotic behavior of the Green function. 
The expansion of logG in powers oll/k (with G defined by (2.7) and (2.5)) is given 
by (6.1). The coefficients of the expansion (equations (6.2) and (4.3)) are expressed in 
terms of the coefficients c n for the expansion of the generalized reflection coefficient, 
which are calculated by using the formula (3.8) with (3.2). The remainder term A N is 
also expressed in terms of c n (see (6.3), (5.3), (4.12), (4.9), (4.6) and (3.9)). Sufficient 
conditions for the validity of the expansion to order N (equation (6.4)) are given by (i), 
(ii), (iii) of section 7. These are not necessary conditions. Equation (6.4) holds under 
broader conditions as shown in section 8. 

We assumed that the potential V(x) is monotone for sufficiently large \x\, but this 
is not an essential restriction for our formalism. The formulas for the coefficients of the 
expansion and the remainder term are valid for any potential, as long as these quantities 
make sense. The particular shape of the potential is relevant only to the conditions for 
the validity of (6.4). The above mentioned assumption on the potential is used only in 
deriving the conditions for (3.12) quoted at the end of section 3. Even for other kinds of 
potentials, we can use the same method to derive the criterion for (6.4). For example, 
in this paper we excluded the cases where V(x) oscillates indefinitely as \x\ — > oo, but 
our method can be applied to these cases as well, if only we study the validity of (3.12) 
(and its analogue for R t ) for such potentials in a similar way as in [19]. 



Appendix A. Properties of scattering coefficients for finite intervals 

Here we summarize some properties of the transmission and reflection coefficients for 
finite intervals. (For details, see [22] and references therein.) First, it is obvious that 



r(x,x;k) — l, R r (x, x; k) — 0, Ri(x, x; k) — 0. (A.l) 
Let us assume that / is piecewise smooth. We have, for finite x and y, 

r(x, y; k) = e lk{x - y) [1 + 0(l/|Jfe|)] , (A.2a) 

R r (x,y;k)=0(l/\k\), Ri(x,y;k) = 0{l/\k\), (A.26) 



as \k\ — > oo (0 < arg/c < n). The first equation of (A. 2b) is a special case of (3.3) with 
N — 1. For finite x and y, equations (A. 2) can be derived more directly from integral 
representations of the scattering coefficients. (See equations (1.14) of [22]. Alternatively, 
we can use equations (3.8) of [22] and the asymptotic forms of the functions and f3 ± 
to derive (A. 2).) 

From equations (3.5) and (3.8) of [22], we can derive the differential equations 
d 

— t(x, y\ k) = ikr(x, y; k) - f(x)R r (x, y; k)r(x, y\ k), (A. 3 a) 

d 

— t(x, y; k) = -ikr(x, y\ k) - f(y)R l (x, y\ k)r(x, y\ k), (A.36) 
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dx 



R r (x, y; k) = 2ikR r (x, y; k) + f(x) 1 — R r .(x, y; k) 



d_ 

dy 

d_ 

dx 
d_ 
dy 



-£-Rr(x, y; k) = -f(y)r 2 (x, y; k), 
dy 



Ri(x, y; k) = -f(x)r 2 (x, y; k), 

Ri(x, y; k) = -2ikR t (x, y; k) + f(y) 1 - R 2 (x, y; k) 



From (A. 2) and (A.3) it follows that, as \k\ ->■ oo (0 < arg/c < n), 

+ 2f(z)e Alk{x - z) + o(l), 

f2/(z)e 4!fc(2 - l) + o(l), 



d 


t(x, z) 


2 

= -2ik 


t(x, z) 




dz 


1 + R r (x, z) 


1 + R r (x, z) 


d 


t(z,x) 


2 

= 2ik 




t(z,x) 


2 


dz 


1 + R l (z,x)_ 


1 + R l (z,x)_ 





t(z,c) 



1 + Rr(z, c) 



T(C,Z 



1 + Ri(c,z) 



dz = —r 

2%k 



dz = — - 

2ik 



r(6,c) 


2 


r(a,c) 


') 


1 + R r (b,c) 




1 + R r (a, c) 





r(c,b) 



1 + Ri(c, b) 



+ 



r c, a 



1 + Ri(c, a) 
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(A.3c) 

(A.3d) 

(A.3e) 

(A.3/) 

(A.4a) 
(A.46) 

(A.5a) 
(A.56) 



Appendix B. Derivation of the short-time expansion 



Let X = x — y. When (9.1) is substituted into (9.4), the integral of the mth-order term 
is proportional to 



oo e i*-X -iwt 



oo K 



m+1 



e -X 2 /(4t) ., 

da; = — 2i F — / 



2t \ m / 2^ \ m ^ 



(B.l) 



where we have changed the variable of integration from u to p = (n/y/t) — [iX/(ty/i)], 
and deformed the contour of integration. The right-hand side of (B.l) can be expanded 
in powers of t by using the formula of Taylor expansion 



(1 + a) 



, m f (m + j-1)! ■ ■ (M + m)! (-l)^c^m 

-^ oJ !(m-l)! 1 j 



(B.2) 



,„.,.,... (M + l)!(m- 1)! (l + 0a) M+m+1 ' 

where M is an arbitrary positive integer, and < 9 < 1. Applying (B.2) to (B.l), and 
carrying out the integration of each term, we have 

xy(4t) ( m_ 2 ™+j(2j + m -l)\ t j+m 



o e iKX e ~iu)t 



oo K 



m+1 



du 



-2i- 



Vt l^? (m-l)!(2i)! (iX)^r. 



+ 



2 m+2M+2 (2M + m + l)\ t M+m+1 /■« p 2M+2 e-P 2 dp 
(m-l)!(2M + 2)! (;x) 2M + m + 2 



/ 



2j/t 



0p 



2M+m+2 



(B.3) 
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From (9.4), (9.1) and (B.3), we obtain the expansion (9.6) with (9.7). The remainder 



where < 9 m < 1. The second term of (B.4) is obviously 0(t N+1 ) as t — > 0. Assuming 
that 1^0, the first term is o(t N ) as t — > if A' N (x,y; k) = o(l/\k\ N ) as \k\ — > oo in 
the region e < arg k < it — e. 
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term A N can be written as 




(B.4) 



